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We consider the generalization of the Kramers escape over a barrier problem to the case 
of a long chain molecule. It involves the motion chain molecule of N segments across a 
region where the free energy per segment is higher, so that it has to cross a barrier. We use 
the Rouse model and find that the free energy of activation has a square root dependence 
on the temperature T, leading to a non-Arrhenius form for the rate. We also show that 
there is a special time dependent solution of the model, which corresponds to a kink in the 
chain, confined to the region of the barrier. The polymer goes from one side to the other by 
the motion of the kink in the reverse direction. If there is no free energy difference between 
the two sides of the barrier, then the kink moves by diffusion and the time of crossing 
t cross ~ N 2 /T 3/2 . If there is a free energy difference, then the kink moves with a non-zero 
velocity from the lower free energy side to the other, leading to t croS s ~ N/y/T. 
The problem of thermally activated escape of a particle over a barrier (the Kramers problem) is very well 
studied (see the review [jl[). Here we consider a generalization of the problem to the motion of a long chain 
molecule over a region where its segments have a higher free energy. The simplest such situation is one in 
which a polymer goes across the interface of two immiscible liquids ||. Recent simulations on this problem 
found that a peptide, when placed in the aqueous phase near a water-hexane interface rapidly translocates 
toward the hexane phase. A more complex and interesting situation is one in which a biomolecule moves 
from one side of a biological membrane to the other jj| . The questions that we study are: given the functional 
form of the barrier, (i) what is the free energy of activation for polymer entry into the barrier? (ii) once a 
part of the polymer has entered the barrier, how much time would it take for the entire polymer to move 
across it? We denote this time by t cross . Muthukumar and Baumgartner Q studied the effect of entropic 
barriers on dynamics of polymers. They found the total passage time to depend exponentially on N . Park 
and Sung || have considered the effect of a small free energy difference between the two sides of the barrier. 
For a flexible molecule they argue that t cross scales as A^ 3 , but this can go over to ~ A^ 2 if there is a free 
energy difference between the two sides. We consider a polymer undergoing activated crossing over a barrier 
whose width w is larger than the Kuhn length I of the polymer, but small in comparison with the length 
of the polymer. Thus, we assume I « w « Nl . In comparison, the entropic barrier Q || has a rather 
large width (~ Nl). We show that for a symmetric barrier, t cross ~ N 2 /T 3 / 2 while if there is a free energy 
difference between the two sides, then t cross ~ N/yT. 

We consider a one dimensional model for both barrier and the chain molecule, described by the continuum 
version of the Rouse model || 

(d t R{n, t) = md nn R{n, t) - V'(R(n, t)) + f(n, t) (1) 

R denotes the position of the n th unit of the polymer at the time t. The boundary conditions for the above 
equation are {9 n R(n,t)} n=0 — {d n R(n,t)} n=N = 0. V(R) is the free energy per unit of the polymer at the 
position R. The portions of the chain over the barrier have a higher free energy per unit. The barrier is 
assumed to be localized near the origin, and to extend in space from — ao to ax, with ao < a\. We take V(R) 
to be a smooth continuous function, with V'(R) — 2kR(R + a ) (R — ai) for — a < R < aiand V'(R) = 
outside these limits. Integrating this and choosing the potential to be zero for R < — ao, we get 

k 

V(R) = -(R + a ) 2 (3R 2 - 2Ra a - ARa x + a 2 Q + 2a ai) (2) 

for — ao < R < a\ and V(R) = -kk (ao — ai) (ao + ai) 3 for R > oi. The barrier height for the forward 
crossing is Vf = jrkao 3 (ao + 2 ax) and for the reverse process, it is V& = ^ka\ 3 (2ao + ai). On crossing 
the barrier, a unit of the polymer lowers its free energy by AV = (ao — a%) (ao + ai) 3 . The form of 
the potential is shown in the figure. Though our results are for this particular form of the potential, the 
conclusions are general and independent of the form of the potential. 
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f(n,t) are random forces acting on the n th segment and have the correlation function 
(f(n,t)f(ni,ti)) = 2^ksTS(n — m)8(t — t\). At equilibrium, the probability distribution functional is 



exp 



^ / dn{±m(<9„i?) 2 + V(R(n))\ 



The free energy of activation for the crossing can be found by 

extremising this functional, subject to the condition that one end of the polymer has crossed R maxi the 
maximum of the barrier. Assuming the polymer to be very long, we find this by taking R{— oo) = — ao 
and the other end of the polymer to be at a point with R > i? max - This extremum configuration sat- 
isfies the Newton's equation md nn R = V'(R), subject to these conditions. For this Newton's equation, 
the conserved energy is E c = ^m(d n R) 2 — V(R(n)) and for the extremum path its value is equal to 
zero. The particle starts at R(— oo) = — ao and ends up at Rf where Rf(> i?max), is the point such 
that V(Rf) = 0. We find Rf = a (-y - y/j 2 - 7) where 7 = (1 + 2^)-. The free energy of activa- 
tion is E act = J R q () ^2mV (R)dR. Interestingly, in the Rouse model, the parameter m is proportional 
to the temperature ||. So if the V(R) is temperature independent, E act a\/T. One evaluation, we find 

77! \fraha\ 



(3 7 2 + 1)VTT3^ - 3 7 (7 2 - 1) In (^7(7 - 1)/ (l + 7 - V^T^t)) 
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The Boltzmann factor e k s T for the crossing of one end of the polymer over the barrier thus has the form 
e -cxms tan t/VT anc j j g independent of N for large N. 

Now, to calculate t cr0 ss, we first look at the mathematical solutions of the deterministic equation, in which 
we replace random noise term f(n,t) by its average. The result is 

(d t R(n, t) - md nn R{n, t) + V'{R) = (3) 

with {d n R(n, i)}„_ = {d n R(n, t)} n=N = 0. The simplest solutions of this equation are: R(n,t) = Rq ( Rq 
is a constant), with Rq < — ao, or with Rq > a\. These correspond to the polymer being on either side of the 
barrier. Thermal noise makes R(n, t) fluctuate about the mean position Rq and this may be analyzed using 
the normal co-ordinates for fluctuations about this mean position. Each normal mode (phonon) obeys a 
Langevin equation similar to that for a harmonic oscillator, executing Brownian motion. Further, the mean 
position (center of mass) Rq itself executes Brownian motion Q . In addition to these two time independent 
solutions, the above equation has a time dependent solution (a kink) too, which corresponds to the polymer 
crossing the barrier. We analyze the dynamics of the chain, with the kink in it, using the normal modes for 
fluctuations about this kink configuration. Our analysis makes use of the techniques that have been used to 
study the diffusion of solitons j?J 

The kink solution may be found using the ansatz R(n,t) = R s {t) where r = n — vt M, when the equation 
(|J) reduces to m d d R £ + vC,^f- = V'{R S ). If one imagines t as time, then this is a simple Newtonian 
equation for the motion of particle of mass m, moving in the upside down potential —V(R) , and v £/m is 
the coefficient of friction. We can easily find a solution of this equation, obeying the conditions R s (t) = — a 
for r — ► -00 and R s (t) = difor r — ► 00. It is R s (t) = (-a Q + e TUJ ( - a ° +ai K 1 ) (l + e™ ( a «+ a ^) _1 , with 
uj = y/k/m and the velocity v — (ao — ai). The solution is a kink, occurring in the portion of the chain 
inside the barrier. The point with r = shall be referred to as the center of the kink. (Actually one has a 
one-parameter family of solutions of the form R s (t + To), where To is any arbitrary contant). As r = n — vt, 
the center of the kink moves with a constant velocity v, which depends on the barrier heights. Note that if 
ao < ai, then Vf < V&, and this velocity is negative. This implies that the kink is moving in the negative 
direction, which corresponds to the chain moving in the positive direction. That is, the chain moves to the 
lower free energy region, with this velocity. If the barrier is symmetric ao = oi( Vf = Vb) then the velocity 
of the kink is zero. 

The center of the kink can be anywhere on the chain and hence the kink can move on the chain. It 
would also execute Brownian motion, due to thermal fluctuations of the medium. To analyze this, following 
'Instanton methods' of field theory jsj, we write 

00 

R(n,t)=R 8 (n-a(t))+^Tx p (t)<l) p (n-a(t),t) (4) 

P =i 

We have allowed for the motion of the kink by allowing the kink center to be at a(t), where a(t) is a 
random function of time which is to be determined. <j) p are a set of functions defined below and X p (t) are 
expansion coefficients. cf> p eventually will turn out to describe motion along the p th normal mode in presence 
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of the kink in the chain. Using the equation (jij) in the equation (|l|) and expanding around the kink, retaining 
first order terms in X p (i), and changing over to the new variables (n, t) with n = n — a(i), (which is just 
segment position along the chain with respect to the center of the kink) we get 

oo oo 

C(u- a (t))dnR s (n) + CJ2 x p (*)4p(M + ^X p (t)C<t> p {n,t) = f{n + a(t),t), (5) 
p— 1 p— 1 

where C = Qd t — mdnn— (a (t)dn + V (R s (nj). We take <f) p to obey the equation 

|Cft — mdnn — (vdn + V (i2 s (n)) j (t) p (n, t) — 0, subject to the boundary conditions cL0 p (n, t) = at both 

ends of the chain, (i.e. at n = — a(t) and at n — N — a(t)). For a very long chain, with the kink located well 
inside it, the ends cannot influence the dynamics of barrier crossing. Therefore, we would make only a minor 
error by imposing these boundary conditions at n = ±N/2. Then, putting p (ri, t) = ^ p (n)e~ Xpt ^~ vi * n l( 2m \ 

we find that ip p (n) has to obey Hip p (n) = \ p il) p (n) with H = |— mfei V (R s (n)) + (t>C) 2 /(4m) j. The 
translation mode of the kink leads to a solution ij>o(n) of the above equation with Ao = 0], Its functional 
form may be found by putting e v< ^/<- 2m ^dnRs(n) =u (a + cuf e^™" (3«o+«i) (i + (a +oi) ) ~ 2 _ Cip (ri), 
where is C the is the normalization factor, chosen such that (ipo \ ipo) = 1. On evaluation, we find C 2 = 
(d^R s (n) je"?™/™! d^R s (n)) = |-vtcj csc(2 tt ) (ai — ao) ao a\. The eigenfunctions ip p (n),p = O..00 form 
an orthonormal set. 

The equation (^) thus becomes 



P =i 



X P (t^pi^e-^^ + (v- a (*)) y(t) = I e < n /( 2m )/(n + a(f), t) (6) 



where Y(t) = ip (n)C + J2p*Li X p {t)e- Xpt ^d^p{n), with <9_ = e^"^ 2 " 1 )^"^™^ 2 " 1 ). On taking the inner 
product of this with tpo(n), we get 



(v-a(t)) lc + J2Xp(t)e~ Xpt/C (i>o \ d^ p {nj) J = (7) 



with £ (t) = i/^ 2 /2 dnVo(n)e^"/( 2m )/(n + a(t),t). Obviously, = and (£o(*)£o(*i)) = <*(* - 

t 1 )(2fc s T/C)/ „ ert , iefcmfc a^e^/™^o(n)] 2 = 5(t-t 1 )fc fl T/(2Cao ax) sec(27rf^) (3 ai -a ) (3a - ai ). 
On taking the inner product of the equation (J^) with ipi(n), we get 



Xi e 



OO 

Xlt/C + (v~ a (t)) J2 X p(t)e- Xpt/C (V* |&J V>p) = &(*) (8) 
p=i 



where £p(t) = | dnip* (n)e v< ' n ^ ( - 2m ' 1 f(n + a(t), t). The equation (||) may be solved, and using the kets 

\X(t)) = X pW I^p) and If (*)> = IV»p>, solution is: 

|X(*)> = \X d {t)) + \X r {t)) (9) 

with \X d (t)) = f [exp (a (f x ) -^A^)}] |X(0)), |X r (t)) 

/* &iT exp { j£ dt 2 (« (*a) - «) 1(*2)} e Stl |£(fi)) ,where A(t) is the operator A(t) = e St ^d-e~ St ^. 
T is the time ordering operator. The equation (^) gives 

a(t) = v + £ (t)/ (C + <^o| fiLe- 5 */^ |X(t))) (10) 

Now taking the average of the above, and using the equation @ and the fact that £o(t) is independent of 



everything that occurred at previous instants, we get ^a (tn — v, so that on an average, the kink moves 
with a velocity v. For the polymer to cross the barrier, the kink has to go in the reverse direction, by a 
distance proportional to the length of the chain. Hence t cross ~ N/v. As v is proportional ^/mk, assuming 
V{R) to be temperature independent we find t cross ~ N/ \fT. 
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If the barrier is symmetric, the kink moves with an average velocity v — 0. If one neglects the kink-phonon 
scattering term in the equation ( |To| ) (i.e., the term (V>o| d-.e~ Ht/,< * \X(t))), then the diffusion coefficient of 
the kink is D = 3|lbT rm (3 ai a ) (3 a - ai ) t , 2 ai^ao \ In the limit t th j f thc 

symmetric barrier as D = \/lF ' ^ n ^ n * s case i the time required for the polymer to cross over the barrier 

is t cross ~ N 2 /D ~ A 2 /T 3 / 2 . 

Park and Sung ^ consider the passage of a polymer through a pore and the barrier is entropic in origin. 
Consequently it is very broad, the width being of the order of N. Hence they consider the movement as 
effectively that of thc center of mass of the polymer which diffuses with a diffusion coefficient proportional 
to 1/A. As the center of mass has to cover a distance A, the time that it takes is proportional to A 3 . If 
there is a free energy difference driving the chain from one side to the other, then the time is proportional 
to TV 2 . In comparison, in our process the barrier is extrinsic in origin and its width is small in comparison 
with the length of the chain. The crossing occurs by the motion of the kink, which is a localized non-linear 
object in the chain whose width is of the same order as that of the barrier. As the kink is a localized object, 
its diffusion coefficient has no A dependence and hence our results are different from those of Park and Sung 
H . In the case where there is no free energy difference, our crossing time is proportional to N 2 (in contrast 
to A 3 of Park and Sung) , while if there is a free energy difference, our crossing time is proportional to A 
(in contrast to A 2 of Park and Sung) . 

It is interesting to note that perhaps, thc kind of behaviour discussed here may have already been observed 
in the simulations of Chipot and Pohoille || . It is also of interest to note that if the barrier height for the 
forward crossing decreased, approaching zero, then the velocity v tends to a finite value. 

I deeply indebted to Professor S.K. Rangarajan, for the encouragement that he has given me over the 
years. I thank Prof. S.Vasudevan for interesting discussions and Prof. B. Cherayil, for his comments on the 
manuscript. 

Figure Caption: The barrier and its inverted form. The barrier heights in the forward 
and backward directions are shown. The dotted line represents the path that determines the 
activation energy 
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